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Abstract 
There appears to be a longtime, very slowly evolving state in dense simple fluids which, 



^ \ for high enough density, approaches a glassy nonergodic state. The nature of the nonergodic 

O 
1^ ■ state can be characterized by the associated static equilibrium state. In particular, systems 



driven by Smoluchowski or Newtonian dynamics share the same static equilibrium and 
nonergodic states. That these systems share the same nonergodic states is a highly nontrivial 
statement and requires establishing a number of results. In the high-density regime one finds 
fi : ... a„ e,.iUb.a«n. =..e. decay, v.a a ..ee-.ep p.oce. ide^iSed . .ode-coupU,, 

theory (MCT). For densities greater than a critical density, 77*, one has time-power-law 
decay with exponents a and h. A key ingredient in identifying the universal nature of 
the long-time results is the recognition that the theory can be expressed in terms of two 
fundamental fields: the particle density p and the response field B. There are sets of linear 
fluctuation dissipation relations (FDRs) which connect the cumulants of these two fields. 
The form of the FDRs is the same for both Smoluchowski or Newtonian dynamics. While 
O ■ we show this universality of nonergodic states within perturbation theory, we expect it to 

m ■ 

be true more generally. 

The nature of the approach to the nonergodic state has been suggested by MCT. It has 

X ■ 
Vh ■ been a point of contention that MCT is a phenomenological theory and not a systematic 

theory with prospects for improvement. Recently a systematic theory has been developed. 

It naturally allows one to calculate self-consistently density cumulants in a perturbation 

expansion in a pseudo-potential. At leading order one obtains a kinetic kernel quadratic in 

the density. This is a "one-loop" theory like MCT. At this one-loop level one finds vertex 

corrections which depend on the three-point equilibrium cumulants. Here we assume these 

vertex-corrections can be ignored. We focus on the higher- order loops here. We show that 

one can sum up all of the loop contributions. The higher-order loops do not change the 



nonergodic state parameters substantially. 

I. INTRODUCTION 

There is a large body of evidence^ that dense simple fluids show an interesting long-time 
slow dynamics. We develop here the ingredients of a theory of the universal features of 
this long-time dynamics. The theory is applicable to systems evolving under Newtonian^"^, 
Smoluchowski^»^ and, possibly, Fokker-Planck dynamics^ in equilibrium. Each of these sys- 
tems share the same thermodynamics and equilibrium static structure. We conjecture that 
there is a long-time slow kinetics regime governed by the solutions of what we call the 
ergodic- nonergodic problem (ENEP). 

In this paper we formulate the ENEP within the self-consistent perturbation theory intro- 
duced in FTSPD^ SDENE^, NDI^, NDII-^, MMSI-^ and SM^. We show, within perturbation 
theory, that both Smoluchowski dynamics (SD) and Newtonian dynamics (ND) lead to the 
same ENEP. We explore the degree to which the long-time kinetics of classes of simple fluids 
depend only on the nature of the interaction potential and not on the details of the kinetic 
process. 

Given the form of an ENEP one must still, for each interaction potential, determine if the 
system undergoes an ENE transition and shows the famous two-step kinetic process explored 
in conventional mode-coupling theory (MCT)^""— . Here we establish the form of the ergodic- 
nonergodic problem for simple fluids. In previous work in SDENE we established the ENEP 
for SD for the case where the dynamical vertices were evaluated using the noninteracting 
ideal gas dynamic three-point vertices. Here we flnd the appropriate ENEP for both SD 
and ND including vertex corrections and determine those properties which determine the 
members of a universality class. 

We show that the ENEP has two sets of contributions. First there is a backbone loop 
structure which includes the one-loop contribution familiar from MCT. Upon this loop struc- 
ture rests a set of vertex corrections. We show here how to add up all of the loop contribu- 
tions. Explicit evaluation of the vertex corrections will be carried out elsewhere. Within our 
systematic development here, MCT is a one-loop theory without vertex corrections. Note 
the vertices in MCT differ from those found here. 

We also show how to practically solve the ENEP in the rather sophisticated approximation 



where one treats all loops but ignores vertex corrections. We call this the no- vertex correction 
(NVC) approximation. In this approximation we removed the constraint that the pseudo- 
potential, used in perturbation theory, be small. 

In working out the NVC approximation in some completeness we have learned about the 
shared kinetic structure that leads to the universal nature of the ENEP. It is interesting 
(but not surprising) that with the proper identification of core variables, (the density p and 
a response field -B), the matrix of cumulants among the variables p and -B, share the same 
fiuctuation-dissipation relations. (The first treatment of these higher cumulants is given in 
NDII.) This is despite the fact that the microscopic definitions of B are quite different for 
SD and ND. What is much less obvious is that the three-point cumulants and vertices for 
both SD and ND obey^ eight new fluctuation- dissipation relations or identities among them. 
These identities are of the same form for both sets of dynamics. We conjecture that one will 
flnd the same set of identities for Fokker-Planck dynamics. This is an exact nonperturbative 
set of results. These statements are true despite the fact the ideal gas realizations of, for 
example, Gb are very different for SD and ND. 

In carrying out the analysis of the self-energies (memory function) and showing that they 
collapse onto a single ENEP, it is crucial that one use these identities and certain theorems 
derived from them. 

Roughly speaking, if one chooses the core variables for each system properly, then cumu- 
lants for two-systems share the same FDRs, thermodynamics, equilibrium static structure, 
and, we show, the same ENEP. 

II. CLASS OF KINETIC SYSTEMS 

We discuss two rather different particle kinetics. For a system of A^ particles with mass 
m, at positions Ri and momenta Pi we can consider Newtonian dynamics, with equations 
of motion 

mki = Pi (1) 

and 

Pi = fi , (2) 



or Smoluchowski dynamics, with equations of motion 

R, = DU + ii . (3) 

The force in both cases is given by 

h = -V.iU , (4) 

where the total potential energy is 

where V{r) is a standard pair potential. D is a diffusion coefficient and ^j is gaussian noise 
with variance 

iUtMt')) = 2DkBT6,J{t - t') . (6) 

Starting with these equations of motion we can reexpress each problem in terms of a path- 
integral formulation. The SD are treated in FTSPD and ND in NDl. Both problems can 
be organised around a field-theoretic treatment of a core problem. All the members of a 
class have the same equilibrium static structure. In particular they share the same static 
equilibrium density correlation functions: 

Cpp,„p{Xi,X2,...,Xn) = {p{Xi)p{x2),...,p{Xn)) ■ (7) 

We now set up the general long-time kinetics in terms of a set of core variables; $ = (p, B). 
The core variables in these systems are the physical observables the number density, 

N 

pil) = Y.5ix,-R,it,)) (8) 

and the crucial response field B{1) which depends on the particle dynamics, SD or ND. 

In FTSPD and NDl we introduced a field theoretical formulation of classical many- 
particle dynamics. The grand canonical partition function is of the general form 

oo N 

Z=Y. W^Tre-^^^-^ (9) 

where the trace is over the phase-space degrees of freedom plus the associated MSR^- con- 
jugate fields, A is the MSR action given by 

^ = ^0 + ^ E "^^^^A- (10) 



where i and j labels space, time and fields p or B, and Aq is the noninteracting MSR action 
and the interacting contribution can be written in terms of the core variables $ which 
essentially defines the response fields B. The interaction matrix aij is given below. Hi is an 
external-time and space-dependent field that couples to the core variables $i. 
The generator of dynamical cumulants is given by 

W[H] = lnZ[H] . (11) 

The cumulants are given by 

55...5W[H] 
^''■■■'~ 6H,6H,...6Hk ■ ^ ^ 

It was shown in FTSPD and NDI that the one-point average 

G^ = ($.) = ^W[H] (13) 

satisfies the fundamental identity 

G, = f r0,e^-<^+^^[^] (14) 

where Tr is the sum over all the defrees of freedom including the conjugate MSR= degrees 
of freedom for a typical particle. For the core variables we write 

N 

where a labels the particles and 

0(°)(l) = 5(xi-i?W(ti)) (16) 

and 0V(1) is strongly dependent on the type of dynamics. For SD dynamics see FTSPD 
while for ND see NDI. Finally we have in Eg. ( 11^ 

AW[H] = W[H + F]-W[H] (17) 

with 

Fi = T.^^J^J (18) 

j 
and where for systems fluctuating in equilibrium the force matrix, appearing in Eq.( !TO|) . is 
given in terms of fourier transforms by 

(^a/siq) = V{q)[SapSi3B + SasS/Sp] ■ (19) 



III. FLUCTUATION DISSIPATION RELATIONS 

A. General Results 

In this section we explore the nature of fluctuation-dissipation relations (FDR) obeyed 
by our class of systems. In ND2 we find the general FDR for Newtonian dynamics for the 
core fields p{x,t) and B{x,t). We found the simple result for n-point cumulants in terms of 
Fourier transforms: 

GBB...Bppil2 ...i,i+l...n)= G;^...^p...p(12 ...£,£+ 1 ... n) (20) 

where 

ipi = B,- ifSuiPi (21) 

and for the single-particle quantities 

gBB...Bpp{i2 ...e,e + i...n) = g;^...^p...pii2 ...e,i + i...n) . (22) 

In the case of SD these relations have been shown to hold for cumulants with up to five 
B insertions^. It seems likely that they hold in general. Do these results also hold for 
Fokker-Planck dynamics? These FDRs are extremely useful. 

B. Two-point cumulants 

Using the FD relations at the two-point level we find for the two-point cumulant: 

GBpil2) = G^,(12) + ^PuJ^G*Jl2) (23) 

where 1 = {qi,uji). Due to translational invariance in space and time 

G^^{12) = Gafs{l)6{l + 2) (24) 



and 



We can also show that 



G'Bp(l) = G^,(l) . (25) 



Gbb(12) = . (26) 



We see that Eq. (l27|) can be written as 

GB,{l) = G*^,{l)+^(5uJlG;^{l) • (27) 



Since Gpp{l) is real, we have the conventional result 



ImGBp{l) = -ImGpBil) = ^^,^(1) (28) 



which can be used to construct the dispersion relation 

du l3uGpp{qi,Q) 



G,b{1) = / — ^ ""^^^'^ . (29) 

^^ ^ J 271 uji- 00 + ir] ^ ^ 

Eg. (12^ leads to the static sum rule 

GpB{qi,0) = -f3Siqi) (30) 

where S{qi) is the static structure factor. Taking the inverse time Fourier transform of 
Eq. ( l29l) we find one of the standard FDRs: 

GpBiq,t) = eit)(3^^Gppiq,t) . (31) 

Combining Eqs.f l27|) and fl25|) we can write the very useful result 

+ i/3uGpp{q, u) = GBpiq, w) - GpB^q, oj) . (32) 

The matrix propagator, Gap, is the first of a number of two-point FDR matrix propaga- 
tors (FDRMP) defined in Appendix B. In addition to Gai3 we have Qai3, Gap, Gap, and 
noninteracting counterparts enter our analysis. 

C. Two-point Vertices 

The two-point vertices Vap.{X) are defined by Dyson's equation: 

Y.Ta,{l)Gpp{l)=5ap . (33) 

This matrix equation reduces, using Eq.( l26|) . to 

rpp(i) = (34) 

TpB{l)GBp{l) = 1 (35) 
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TBpil)G,Bil) = 1 (36) 

and 

rBB{l)GBp{l) + TBp{l)Gppil) = . (37) 

Eq.f l37|) can be written in the form 

G,,{l) = -G,B{l)rBB{l)GBp{l) . (38) 

We easily find 

rBp(g,w) = r*5(g,w) . (39) 



Starting with Eq.f l27j) . we can use Eqs. fH2|) and f l36|) to write 

1 1 



+ tl3iOiGpp{l) . (40) 



T,b{1) Tb,{1) 
Using Eq.f lSSp and canceUing a common denominator, gives 

i^u^YBBil) = r,B{l) - Tbp{1) (41) 

which can be compared with Eq. fl32|) . 



D. FD Relations and Three-point cumulants 

There are substantial differences in specific cumulants for the different realizations of 
simple fluid dynamics. For example the noninteracting three-point cumulant G^^ (123) for 
SD and ND are quite different. 

The FD relations, Eq.(l20l). for three-point cumulants, reduce to eight independent iden- 
tities: 

GbBB = G*BBB + il^^lG*pBB + il^^2G*BpB + il^^zG*BBp (42) 

- /3^ {uJiUJ2G*ppB + ^l^3G*pBp + UJ2UJ-iG*Bp^ - il3^UJiU2UJ^G*ppp 

Gbbp = G*j,sp + iPu^G*pj,p + iPu:2G*j,pp - /^'c^ic^aG^^p (43) 



Gbpb = G*BpB + iP^iG*ppB + il^^zG*Bpp - l3^ujiujzG*ppp (44) 

GpBB = G*pSB + ^P^2G*pps + iP^zG*pSp - fi'^2UJzG*ppp (45) 

Gbpp = G*^pp + iPu,G;pp (46) 

GpBp = G;j,p + i(5u2Glpp (47) 

GppB = G;p^ + ^^u:,G;pp (48) 

Each of these constraints give an identity relating the 3-point cumulants. Since we know 
Gbbb = G*BBB = 0) Eq.dH]) reduces to 

= il3ujiG*pBB + iP^2G*BpB + iP(^zG*BBp 
- 13'^ {uJiUJ2G*ppB + UJiUJ^G*pBp + UJ2UJ^G*Bpf)j - i/3^UJiUJ2UJ3Gppp . 

As shown in ND2, these cumulant identities lead to useful relations satisfied by the three- 
point vertices. The relationships connecting the three-point cumulants and the three-point 
vertex functions are given by 

r,^,(123) = - Y. r„,,(l)r^^,(2)r,,,(3)G,.^v'(123) (49) 

a'fi'i/' 

If we introduce some simplfying notation: 

01 = if3uiTBpp (50) 

02 = i(3u2TpBp (51) 

03 = i/3u3TppB (52) 

V'l = /3^UJ2UJ3TpBB (53) 

^2 = /3^ujiU3Tbpb (54) 
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ips = /3^u;iU2Tbbp (55) 

P = -i/3^UJiUJ2UJ3TBBB , (56) 

then we find, after algebra, tlie FD vertex relations: 

P + ^3 + ^2 + ^1 + 03 + 02 + 01 = (57) 

- 0* + 01 + 02 + ^3 = (58) 

-0; + 0i + 03 + V2 = O (59) 

-0t + 02 + 03 + ^-1 = . (60) 

IV. KINETIC EQUATION 

Now we start working toward a description we can use to identify the ENE problem. As 
a first step we want to obtain the form of Ta/^iq, oo) in the time domain. If we look at 

and note that the response functions vanish algebraically as cj — ?> oo, we see that the two- 
point vertex diverges in this limit. We assume that 

lim Tspiq, w) = -DqCj^ + icuAa + Cq (62) 

where the coefficients A, C and D are real. Define the "local" quantity 

T%{q, to) = -Dqco^ + icoAq + Cg (63) 

for all frequencies. We also define the subtracted quantities: 

rg(g, u) = TbM, CO) - T%iq, u) (64) 

which vanishes for large frquencies. At low frequencies 

T%iq,0)=TBpiq,0)-T%iq,0) (65) 
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where, using Eqs.( l30|l and ([6T 



TB,{q^O) = -{PS{q))-' 



(66) 



and 



r%{q,o) = c,. 



This leads to 



(67) 



<s) 



rS(g,o) = -(^5(g))-i-c, 



We assume that the FDR holds for T^^'^ and 



.{£) 



PuT'^>siq,0) = -2ujA, 



(68) 



(69) 



or 



PT^^Uq,0) = -2A, 



(70) 



and we assume Ag > 0. 



Taking the inverse Fourier transform of Eq. fl63|) gives 



and 



rg(?,t-t') 



— - A ^ 



-^oTTT^T ~ ^aTTT + ^q 



5{t - 1') 



(71) 



.W 



PV%{q,t-t') = -Ag5{t-t') 



(72) 



We have then 



.W 



W, 



rBp(g, t - t') = r^;,(g, t-t') + T%{q, t - t') 

TsBiq, t-t') = -r^Ag5{t - t') + rgjj(g, t-t') . 

It then follows for the subtracted parts that we have the dispersion relation 

du Pu:T%\{q,uj) 



<s) 



n'M^^) 



27r 00 — 00 + ir) 



(73) 
(74) 



(75) 
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reflecting the fact tliat Fi^ is analytic in the upper-half plane and vanishes as a; — t- oo. In 
the time domain the FDR reads 



^W, 



,'\ ^ r(«) 



rS(g, t-t') = m - t')^r'BUq, t-t') ■ 

Now express Eq. fl37|l in the time-domain 

where the g-dependence is surpressed and 

^(t, t') = f dsT.B,{t - s)Gpp{s -t')+ f dsT.BB{t - s)G,B{t' - s) 

J — OO J — OD 

The self-energies are defined by 



(76) 



(77) 



(78) 



(79) 



and we have used the fact that T,Bp{t — s) ^ 9{t — s) and GpB{t' — s) ^ 9{t' — s). We then 
use the fluctuation dissipation relations 



^Bp{t -S)= e{t - s)^l3T.BB{t - S) 



GpBit' -s) = 9{t' - s)^f3Gpp{t' - s 



(80) 



(81) 



to obtain 



^{t,t') 



ds 



d_ 
ds 



l3T.BB{t-s) 



d 



Gpp\s — t ) 



- / dsJ:BB{t-s)—l3Gpp{t'-s) . 

J—OO US 

Integrating by parts in the first integral: 

^{t,t') = ~^j:BB{0)Gpp{t-t')+ f ds/3^BB{t-s)-^Gpp{s-t') 

J—OO OS 



(82) 



d 



I dspi^BBit- s)—Gpp{t' - s) 

J—OO OS 



d 



-l3Y.BB[Q)Gpp{t - t') + j^ ds/3EBB{t - s)—Gpp{t' - s 



(83) 
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where we assume t > t' . We then have, putting in the g-depence, the kinetic equation 
f g!,(g, t)Gpp(g, t-t') = -P^BBiq, t = 0)G,,(g, t - t') 

+ J^^ ds/S^BB {q,t-s)—Gpp{q,t'-s) (84) 

where for t > t' we can drop the TB^{q)GBpiq,t — t') term. Putting all these results together 
we have the kinetic equation 






Gpp{q,t)- J^^dsPJ:BB{q,t-s)—Gpp{q,t'-s) = . (85) 



We see that our dynamical problem is now in the form of a memory function equation 
and the dynamic part of the memory function is given by the self-energy T^bb- The static 
structure factor is the same for all fluids essentially by definition. The coefiients Ag and Dq 
can be constructed using perturbation theory. Let us look at the two cases studied so far. 
At zeroth-order the two-point vertex is given in the SD case by 

(0)/ X i-iu + Dq^) 
lVM.-) = j^^^ ■ (86) 

From Eq.( l63|) we can read-off 

Df^iq) = (87) 

and 

In the ND case we have the more complicated result: 



where 



,2 



and 



S{z) = 1 - 2ze-' / dwe" - i^/^ze-' (91) 

^0 



z = -^ (92) 



V2qVo 
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with mVQ = ksT. In the large frequency hmit 



.(0) 



TBpiQ,^) 



f5-p 



and we can identify 



i)^' 



1 / 1 \' 



"°^^^^^ ^ -^ u 



<-(g) = 



1 



cri^)-j-^ • 



(93) 

(94) 

(95) 
(96) 



The kinetic equation, Eq.( l85l) . is diagonahzed using a Laplace transform in time. We 
define 



and we need the results 



and 



Gf,f,{q,z) = / dte '*Gpp{q,t) 
Jo 



dte ''Gpp{q,t) = zGpp{q,z) - S{q) 



dte-^'Gppiq,t) = z[zGppiq,z) - 5(g)] 



(97) 



(98) 



(99) 



Laplace transforming the kinetic equation gives 

D{q)z[zGpp{q, z) - S{q)] - A{q)[zGpp{q, z) - S{q)] - ^^ S-\q)Gpp{q, z) 

- l3f:BBiq,z)[zGpp{q,z) - S{q)] =0 

Dividing by [zGpp{q,z) - S{q)] gives 

p-'S-\q)Gppiq,z) 



D{q)z-A{q) 



/3i:BB{q,z) = . 



(100) 



(101) 



[zGppiq, z) - S{q)] 

We assume in the long-time limit and for high densities that Gpp{q,z) and T,bb{(1,z) blow 

up as 2; — )■ 0. So the D{q)z — A{q) term can be dropped. The kinetic equation, Eq. OlOll) . 

reduces to 

Gpp{q,z) 



-P'S{q)^BB{q.z) 



(102) 



[zGpp{q, z) - S{q)] 

We must now turn to the mechanism which produces large Gpp{q,z) and TiBB{q,z) as 
z -^ 0. This involves determining E^b as a functional of Gpp. Solution of the resulting 
self-consistent equation, Eq. (ll02p . for Gpp corresponds to what we call the ENE problem. 
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V. DETERMINATION OF SELF-ENERGIES 

A. Dynamic Ornstein-Zernike Relation 

We now want to construct the self-energy T,bb self-consistently as a functional of Gpp. 
The first step is the the derivation^ of a Dynamic Ornstein-Zernike Relation. 
Inserting Eq. (IT^ in 

r - ^ r 

we can use the chain-rule for functional differentiation to obtain 

Gij = Qij + 2_^ CikGkj (103) 

k 

where the single-particle FDRMP is given by 

g,, = rr0,0,e^'^+^'^ (104) 

and 

Q,- = tr^.e^-^+^^-^AW^. (105) 

OLrj 

Since AW can be treated as a functional of Gi we see at this stage that we have available a 
self-consistent theory. If we define the matrix-inverses 

Y.TikGkj = 6,, (106) 

k 

and 

^likGkj = Sij (107) 

k 

then the two-point vertex is given without approximation by 

r^,- = jij + K,, (108) 

where 

Kij = -J2'yikCkj (109) 

k 

is the collective contribution to the two-point vertex. 
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B. General Analysis 

Much of the analysis of the two-point vertex can be carried out to all orders in an 
expansion in the potential. In general the two-point vertex is given by 

J- ij = lij ~ O'ij — l^ij 

where the collective self-energy, Sjj, is given by Eqs. fllOQp . f llOSp . and f ll07p . (sums over 
repeated indices implied) 

"*-^i n=2 ^• 

where we have used Eqs.f ll7p and ( [T2|) and generated the functional power-series in Fj. Since 
Fi is defined by Eq. fl38|) and is independent of Gi, we have 



oo -1 r 

n=2 ^- "^i 

After defining the single-particle quantity 

Gkviv2...v„ = Tr0fce^''^+^'^0„i0^2 ■■■(Pv„ (111) 

and using the chain-rule for functional differentiation, we find 

OO -| 
^ij lik / , 7ykv\V2...Vn'^U-iV\'^U2V2 ■ ■ ■ '^ UnVn^ U\U2 . . ■Unl'- Ij ■ \^^'^) 

n=2 '^■ 

Introduce the amputated vertices V and 7: 

ykviV2...Vn ykmyvtwiyv2W2 ■ ■ • yvnWnimwiW2---Wn 1-'--'^'^/ 

and 

^kviV2...Vn ^km^V\Wl^V2W2 • • • ^VnWn'- mWlW2---Wn ■ \^^^) 

Then we have quite generally 

00 -| 

^ij Jik / J 7ykmyviwiyv2W2 ■ ■ ■ yvn,Wn,imwiW2---Wn 

n=2 ^• 

^ '^u\Vi'^U2V2 • • • ^u„v„'~J^ uixi'~^ U2X2 ■ ■ ■ ^u„x„'~^p£^ pxiX2---x„^ £j ■ \^^'^) 
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We see factors of the FDRMP 

uv 

entering the development. We can then write 

OO -| 
^ij / y 7liwiW2...W„^XlWl^X2W2 ■ ■ ■ ^XnWn'- JXIX2...X„ ■ \^^') 

n=2 ^• 

In terms of Fourier transforms we have the general result: 

OO 1 „ 

S„^^^, (12) = E - / d2d3 ...d{n + 1)5(1 + 2 + . . . + (n + 1)) 

71=2 ^- '' 
X 7a,ai...a„(12 . . . (u + l))^^,^, (l)^^^^^ (2) . . . Ga^^^^^^^„^,^{n + 1) 

xf^^.^,...^„(12...(n + l)) . (118) 

Now we are interested in the BB component 

°° 1 f 
Sbb(12) = y — dld2 . . . dn6(l + 2 + . . . + (n + 1)) 

X 7Bai...a(„+i) (12 . . . (n + l))Ga,^,{l)Ga2^2 (2) • • • G^f^+DMCn+i) (^ + 1) 



X 



fBMi...Mn+i)(12...(n+l)) . (119) 



In the long-time regime the internal propagators are dominated by pp lines and we have in 
the low-frequency regime 

OO -, 

Sbb(1) = E ^ / d2d3 ...d{n+ 1)6(1 + 2 + . . . + (n + 1)) 

X 7V..p(23 ■■■{n + l))G,,i2)G,,{3) . . . G,,{n + l)f Bp...p(23 . . . (n + 1)) . (120) 

Next we note that the major contribution to the integals is from the low-frequency portions 
of the Gpp. In that regime we can evaluate the vertices at zero-frequency. At zero-frequency 
we can use the vertex theorems (see Appendix A) 

7Bp...p(gi, 0, ^2, . . . kn+i, 0) = --— (121) 

pp 



to obtain 
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f Bp...p(gi, 0, ^2, . . . K+i, 0) = -/3-^7(gi, /ca, . . . , k^+i) (122) 

/3'Sbb(1) = E ^ / ^2d3 . . . d(r2 + 1)5(1 + 2 + . . . + (n + 1)) 

X 1(5,^(2)^^^(3) . . . Gpp(n + l)7(gi, ^2 . . . A;„+i) . (123) 

This is the long-time approximation for the collective part of the memory function. If we 
ignore vertex corrections and replace the static vertices with its ideal-gas form 

7(giA:2 . . . A;n+i) = — (124) 

pn 

we obtain 

(5^T.bb{1) = E ^ / ^2rf3 ...d{n+ 1)5(1 + 2 + . . . + (n + 1)) 
n=2 '^^ -^ 

xi^G,,(2)G,,(3)...G,> + l) . (125) 

If we express the Gpp{ki,Ui) in terms of its space-time Fourier transform we can do the 
frequency integrals and obtain the very simple result 

OO 1 

/32Sbb(1) = Y.-J dh^^'"^' j d''ne''--''^{-D{nM)T (126) 



n=2 

where we have 



Gppijri^tij 



D{nM) = '^1^^^^^ . (127) 

We can do the "loop" sum to obtain 

/32Sbb(1) = / dUe''^^^ J d^ie'^'-'^HlDin, t,)] (128) 

where 

H[x] = e-'' -1 + x . (129) 

If we expand in powers of D and keep the lowest nonzero result we obtain our previous 
perturbative results at second order—. We are able to sum up all the loop diagrams, wihout 
vertex corrections, and go well beyond MCT. 
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We also assume, consistent with the long-time approximation, that for low-frequencies 

G,,{q,u):^G,Biq,uj) (130) 

and in Gpp{q, u) the Gpp{q, u) term dominates and its coefficient can be replaced by its w = 
value. Then using the FDR we find 

^^^^ = -/3pV{q)Gpp{q,u)/p' 

= -V{q)S{q)F{q, u)/p = -D(g, u) (131) 

where 

Viq) = (3pV{q) , (132) 

Gppiq,oj) = Siq)Fiq,oj) , (133) 

and 

S{q) = pS{q) . (134) 

Taking the inverse Fourier transforms we obtain the very simple result 

/3^j:BB{r,t) = H[D{r,t)] . (135) 

Finally, in conventional MCT the ENE problem is formulated entirely in terms of the 
structure factor. In our formulation here this is true only if we ignore vertex corrections. In 
general one must also include the n-point vertex corrections at n-loop order. This will be 
persued elsewhere. 

VI. ANALYTIC LONG-TIME ANALYSIS 

Putting the Fourier-Laplace transform of Eg. (11351) into Eq. fll02p . we find that the univer- 
sal long-time kinetics is governed in the no- vertex correction approximation by the equation: 

Fiq,z) 



l-zF{q,z) 



//■OO 
£re''^-' dte-''H[D{r,t)] (136) 
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where 

H[x]=e-'' -1 + x (137) 

and 

D{r, t) = -_J ^,e-^'"V{q)S{q)F{q, t) . (138) 



Fourier-transforming over space 



D{q,t) = -V{q)S{q)F{q,t) (139) 



or Laplace transforming over time 



D{r, z)=^-j ^e-'^-\/(g)5(g)F(g, z) (140) 

or the double transform 

D{q,z) = -V{q)S{q)F{q,z) . (141) 

If we replace H[x\ with ^x^ we have the ergodic-nonergodic problem studied in SDENE, 
ND2 and SM. Using the full H[x\ we have the ENEP where we keep all loops. Note the 
ENEP is the same for SD and ND. 

Now we look for a solution of Eq. (ll36p of the form 

F(r,t) = /(r)+V'(r,t) (142) 

or 

F{q.t) = f{q)+ij{q,t) (143) 

valid at long times and where ip{r^ t) is "small". f{q) is the crucial nonergodicity parameter. 
In terms of the Laplace-transform over time 

fir) 
F{r, z) = ^^-^^ + ij{r, z) (144) 

and the Fourier transform over space 

F(g,^) = ^+^(g,z) . (145) 
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The quantity D{r,t) can be written 



D{r, t) = -_J -^,e-^'-'V{q)S{q) [f{q) + ^(g, t)] 



where 



and 



L>o(r) + A(r,t) (146) 



^o(r) = -_j ^,e-^'-V{q)S{q)f{q) (147) 



A(r, t) = 1 1 ^e-«-V(g)S(g)V^(g, t) (148) 



or 



A(r, z) = - / -^e-^'''-F(g)^(g)^(g, ^) . (149) 

Expanding in powers of A we find keeping terms of second order 

/f[D(r,t)] = iJ[Do(r) + A(r,t)] 

= i7[Do(r)] + A(r,t)(l-e-^«('")) + -A2(r,t)e-^o(^) . (150) 

This is inserted into the right-hand-side of Eq. 0136p . The left-hand side of the same equation 
takes the expanded form 

F{q,z) 



L{q,z) 



l~zF{q,z) 



f{q) ^ ip{q,z) ^ zij^{q,z) 



z{l-f{q)) {l-f{qW (1 -/(?))' 
Eq. (ll36p takes the form 

L(g, z) = R,{q, z) + R^^q, z) + R^{q, z) (152) 

where 

R,[q^z) = S{q)^^ (153) 
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with 



and 



Ho{q)= f d^re'''-'H[Do{r)] (154) 

i?2(g, z) = Siq) I £re'''-'A{r, z){l - e'^''^'')) (155) 

//"OO 1 

dVe*'?-'" / dte-"*-A2(r,t)e-^«('^) . (156) 

JO 2 

Now we match orders in our perturbation theory. The leading order goes as 1/z for small 

z. Matching coefficients at leading order we have 

m 



and 



1 - fiq) 

and the higher-order contributions are given by 



S{q)Ho{q) (157) 



ij{q,z) z^^{q,z) D ^ N , D / A /1e;s^ 

+ 71 — 7rT77 = ^2(g,^) + ^3(g,^) • (158) 



Eq. (ll57p is critical in the development. The question whether we have an nonergodic phase 
is determinecl= by this " static" quantity. A nonzero / depends on static parameters. 

Let us focus on the higher-order contributions to Eq. (ll58p . Consider the hnear contribu- 
tion. Using Eq. (ll55p and Eq. (ll49p . we have 

J^C,,kHk,z) (159) 

where we have defined the important matrix, 

Cg^^ = S{q)- /'dVe*'-''(l - e-^'>^''^)e~'^-^V{k)S{k) 
P '' 

= S{q)M{q-k)V{k)S{k) (160) 

and where we have 

M{q) = /dW'?"(l - e-^°('')). (161) 



The higher-order contributions, Eq. (11580 . can now be written 
ipi^q^z) zip'^{q,z) f (Pk 



+ n - mP - 1 ll^^'-^'^'"- -'' + ^^'■'- '1 ■ ('«') 



(1 - my (1 - /(?))' J (21) 



23 

It is notationally advantageous to introduce 

^{q,z) = il-fiq)f<P{q,z) (163) 

and find in terms of </> 

^—^^C,,,<P{k,z) + R,{q,z) (164) 

where 

C,,k = S{q)Miq~k)Vik)Sik)il-f{k)y . (165) 

Let us turn to Rs^q, z). Starting with Eqs. (11561) and (11481) we have: 

i?3(g, z) = Siq) [ dhe"'-'' r dte-''lA^{r, t)e-^'^''^ (166) 

J JO 2 

and 

1 r d^k 
A{r,t) = -pJj^s^~""'nk)S{k){l - mn{k,t) (16?) 

we can write 

where 

C,,k,P = l^|rfVe*(^-'^-^)-^e-^"MF(A;)5(fc)(l - f{k)fV{p)S{p){l - f{p)f . (169) 

Pulhng all of this together we can rewrite Eq. (ll62p as 

d^k 



{2n) 



:Kk - C.^Mk, z) + z{l - f{q))<P\q, z) 



d^k f d^p ^ /•» _,^ _,, 



(27r)3y (27r)3^^'''^io ^^^ <^(^'^)<^(P'^) • (I'^O) 

Continuing our perturbation theory and matching equal powers, we assume 

<P{q,t) = ^^'\q,t) + <P^'\q,t) (171) 

and require that the first-order terms must satisfy 

J-0^,i^,k-C,,,)<p^'\k,z) = O . (172) 
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Then at second order we have 



{S,,k - C,,,)0(2)(A;,z) + z{l - f{qm^'\q,z)f 



{2tt 

<fk f d^p 



/•LXJ 

Cg,k,p dte~-^'ct>^^\k,t)ct>^^\p,t) . (173) 



{2nfJ {271)^'^"''' Jo 
Looking for a first-order solution of the form 

<P^'\q,t) = Ae.Mt) (174) 

where A is independent of q and t, and e^ is the right eigenfunction of Cq^k with unit 
eigenvalue 

{5g,k - Cg,k)ek = . (175) 



(27r)3 

We also introduce the left eigenfunction 

d^k 



(2 
We normalize the eigenfunctions using 



/7^^('^.,.-^.>';^ = . (176) 



d^k 



(2vr) 
and 

d^k 



/(WJ--' (ITT) 



(27r)3 

The second-order contribution is then given by 

d^k 



e,il-fik))ei = l . (178) 



(27r: 



AS,,, - C,,,)(t>^'\k,z) + z{l - f{q))[Ae,<PAz)] 



2 



d k f d p ^ ^ ^ a2 I j^-zt,2 



oo 



Matrix multiply by tg. This kills the first term on the left and after using the normalizations 
gives the rather simple result: 

roo 

zA^cPliz) = A^\ / dte~''(f)l(t) (180) 

Jo 
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where we introduce the parameter 

<fq r <fk r d?p 

(2^ J (2^ J (2^ 

and Eq. dlSOp reduces to 



r d-'q r d-'k r dy 
^ - i 72^ 7 72^ i 7^^e,C,,,,,e,e, (181) 



z(j)'^{z) = X / dte-''(j);{t) . (182) 

Jo 

Eg. (11821) is satisfied by a power-law solution 

Mt) = t'" (183) 

with Laplace transform 

^"W = ^^ • (184) 

Putting this into Eq. (I182p we find 

or, the MCT result, 

g^- ■ 

As discussed in SM, one can go further and treat the full two-step decay process including 
the von Schweidler— contribution. 

VII. PSEUDO-POTENTIAL 

In SDENE we worked to second order in the pseudo-potential V{q) which was expressed 
in terms of the structure factor S{q). The development was self-consistent and thorough. It 
involved the construction of the second-order self energies S^ ' and S^^ , showing that 
these satisfy the appropriate FDR and the associated static sum-rule giving the coorespond- 
ing static approximation in SDENE. 

In the present case it is easier to simply work out the static manipulations giving the 
direct-correlation functions, Cniq), as a power-series in the static pseudo-potential. In the 
all-loops approximation we ignore vertex corrections and do the loop sums to obtain the 
static approximation 

- pCoiq) = V(q) _ ^ /"^ rdr sm(qr) H (Y (r)) (187) 

q Jo 
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I POO _ _ 

Y{r) = — / qdqsm{qr)V{q)S{q) (188) 

VlTirjr Jo 

and 

H{x)=e-'' -l + x . (189) 

VIII. ENE PROBLEM 

One then has the ENE problem to solve. Given a static structure factor, is there a 
solution to Eg. (11571) with /(g) ^ 0. This requires a numerical treatment of Eg. (11571) . For 
those cases where we find a nonergodic regime, f{q) ^ 0, we must solve the eigenvalue 
problems posed by Eg. (II 75 p and Eg. (II 76 p with the matrix Cq^k defined by Eg. (11650 . Then 
one is in a position to compute the parameter A given by Egs. (11801) and (11690 . The power-law 
time-decay exponent is given the by Eg. (1 186 p . This numerical analysis has been carried out 
in the simplest case of hard-spheres and using the Percus-Yevic approximation for the static 
structure factor. An ENE transition is found to occur for reduced density rj* = 0.651.... A 
full numerical analysis will be given in a separate paper. 

IX. SUMMARY 

We have shown here how one can the study of the ENE transition for dense fluids well 
beyond the one-loop treatment of mode coupling theory. We have shown how one can treat 
multiple-loop calculations in a straight-forward manner. We focussed on the approximation 
where we ignore static vertex corrections and sum all loop contributions. Still left to be 
worked out is the role of vertex corrections. Do they change the picture developed here 
in any significant way? This can be tested by doing perturbation theory in the pseudo- 
potential. This will be discussed elsewhere. 

A key result is that we have shown that SD and ND share the same ENE problem. It 
would be interesting to investigate whether the more elaborate Fokker-Planck dynamics falls 
into the same universality class. We speculate that it would lead to the same ENE problem. 
As a first step one would need to show that the core-variables p and B satisfy the same 
FDRs as for SD and ND. 
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Appendices 

X. HIGHER-ORDER THERMODYNAMIC SUM RULES 

A. B-Theorem 

In this section we show how the fluctuation-dissipation relations (FDR) can be used to 
derive thermodynamic sum rules. We start with the general FDR for the core fields p(x, t) 
and B{x,t). The n-point cumulants in terms of time- Fourier transforms obey 

Gbb...Bp...p(12 ...i,i + l...n) = G;^...^,...,(12 . ..£,£ + 1 .. .n) (190) 

where 

ipi = Bi- if3uiPi . (191) 

Let us now consider the case of n 5- variables and one p variable. We can write 

G'pB...B(12...n + l) = G;^ .^(12...n + l) . (192) 

Let us introduce a little notation. If we remember the conserving frequency ^-function we 
can write 

GpB...B{l'2 ...n + l) = GpB...B{u;i,U2, ..., a;„, Un+i)S{ui + U2 + ■ . . + ojn+i) (193) 

where we surpress the wavenumber dependence. Using this notation in Eq. (I190p and setting 
U2 = oj-s = . . . = Un = 0, we have cj^+i = —Ui, and 

GpB...B{^i, 0, . . . , 0, -wi) = G;^, j,{uji, 0, . . . , 0, -ui) - il3uJiG*pB^^p{uJi, 0, . . . , 0, -a;i)(194) 

or 

2^G';^...b(^i, 0, . . . , 0, -u{) = -z/3a;iG';5...^(u;i, 0, . . . , 0, -uj{) . (195) 

Notice that 

GpB...Bp(^i,0,...,0,-a;i) = . (196) 

Now consider the Fourier representation: 

(^pB...B{ti,t2,...,tn+l)-J-r—e J ~2~ -*■'■■■ 
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X / ^e-^-"(*-*"+^)G,B...B(a;i,a;2, • • • ,a;n+i) (197) 

where oOn+i = -^i - CJ2 - • • • - Wn- Now if t„+i > ti then (^^^...^(ti, t2, . • . , t„+i) = 
0. This imphes that Gpb...b{^i,^2, ■ ■ ■ ,^n+i) is analytic in the upper-half plane. So 
GpB...BB{^i, 0, . . . , 0, — wi) is analytic in the upper half plane. We can then construct the 
dispersion relation 

GpB...BB{(^i,0,---,0,-uji} = ^^ : . (198) 

J 71 UJ — UJi — 17] 

We showed earlier 

G"pb...b{^i. 0, . . . , 0, -c^i) = -^G;b..M^ 0, . . . , 0, -c^i) (199) 

so we have 

GpB...BB{'^l,0,...,0,-Ui) = - ^-— ^^ : . (200) 

J 2n u — Ui — IT] 

Taking cui = we find 

GpB...BBiO, 0, . . . , 0, 0) = - 1 ^f3G;j,,,,spi^, 0,...,0,-u) . (201) 

Returning to the time domain we have 

^pB...p\tl,t2,- ■ ■ .tn+l) — — — e I ~0 — -^ ' ■ ■ ■ 

X / ^e-^-"(*"-*"+^)G,B...p(c^i, c^2, . . . , c^n+i) . (202) 

Fourier transforming over t2, ■ ■ ■ ,tn we find 

GpB...Bp{tl,(^2, ■ ■ ■ ^n, tn+l) = 



-e 



i.;i(ti-t„+i)e^Wi(-^2+...+^n)^^^^^^^^(^^^0,...,0,-Wi) . (203) 
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Taking all the external frequencies to zero we obtain 

dui 



GpB...Bp{tu 0, . . . , 0, tn+l) = J ^e-*-^(*^-*"+i)GpB...Bp(c^i, 0, . . . , 0, -Ui) . (204) 



Setting tn+i = ti we obtain the important result 

dui ~ 



GpB...Bp{ti,0,...,0,ti) = J^GpB...Bp{ooi,0,...,0,-ui) . (205) 
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Using this result in Eq.f l20ip we find 

G'pB...BB(0,0,...,0,0) = -/3G,B...Bp(ti,0,...,0,ti) . (206) 

If we rewrite tliis in tlie form 

GpB...BB(0,0,...,0,0) = -/3G'ppB„,B(ti,ti,0,...,0) . (207) 

It is easy to see that one can turn the crank and find 

GppB...Biti,ti,0,...,0) = —f3GpppB...B(ti,ti,ti,0,...,0) . (208) 

and 

GpB...BB{0, 0, . . . , 0, 0) = (— /3)"'G'pp...p(ti, ti, . . . , ti) 

= {-f3rSc{qi,q2,-.-,qn+i) (209) 

where Sc is the n + 1 static density cumulant. 

B. n-point Amputated Vertices 

We define n-point amputated vertices with 

(:?„,,,...„„ (12 ...n) = G^,p,{l)G^,,,{2) ■ ■■G^„,M^,,p,,„p^{12 ...n) . (210) 

Notice that 

Gbb...b(12 ...n) = GBpil)GBpi2) ■ ■ ■ GBpin)Tp,p„„p„ (12 . . . n) = . (211) 

We are interested in 

GpB...B{^2...n) = GpBil)GBpi2)---GBpin)TBp...pil2...n) . (212) 

Setting all external frequencies to zero and using Eq. ( 1209^ cancelling the common factors of 
— /3 we find the very useful result 

- /3rBp...p(0, 0, . . . , 0; gi, g2, • • • , Qn) = lpp...p{qu ?2, • • • , qn) (213) 

where 7 is the static n-point amputated vertex. 
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C. Single-particle Quantities 

It should be noted that this analysis goes through for the single-particle quantities. Thus 
we can replace GpBB — ^ QpBB- Then we have the identities 

Gp{xrM) = Tr0,(xi, ti)e^^ = p (214) 

Qpb{xi,X2]Q) = -f3Tr(f)p{xi,ti)(f)p{x2,ti)e'^^ 

= -(5p5{xi - X2) . (215) 

QpBb{xi,X2, X3; 0, 0, 0) = f3'^Tr(f)p{xi,ti)(f)p{x2, ti)0p(a;3, ti)e^^ 

= /3'^pS{xi — X2)S{xi — xs) (216) 

GpBBBixi,X2, X3, X4; 0, 0, 0, 0) = -/3^Tr0p(xi, ti)0p(a;2, ti)(f)p{x3, ti)(f)p{x4, ti)e'^^ 

= —(3^p5{xi — X2)5{xi — x^)5{xi — X4) (217) 

an so on. Fourier transforming over space: 

QpB{qi,q2\Q) = -Pp5{qi + q2) (218) 

GpBB{qi, q2, qs; 0, 0, 0) = f3^p6{qi + q2 + qs) (219) 

GpBBBiqu g2, gs, qr, 0, 0, 0, 0) = -/3^p6{qi + ga + gs + ^4) (220) 

and so on. Then, for the associated amputated vertices, 

iBpiquO) = -j-_ (221) 

lBpp{qi,q2,q3;0,0,0) = ^ (222) 

iBpppiqi, g2, gs, g4; 0, 0, 0, 0) = -— . (223) 
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XI. FDR MATRIX PROPAGATORS 

A FDR matrix propagator (FDRMP) A^i^(q,u) satisfies the properties 

A^,{q,u) = A:^{q,u) (224) 

+ i(3ujApp{q, u) = ABpiq, u) - ApB{q, u) (225) 

ApB{q,uJ) = / 7^ , . 226 

J 271 U — U + IT] 



and 



ApB{q,0) = J^f3App{q,u) (227) 



with 



ABB{q,oo) = . (228) 

If Aai3{q,uj) and Cai3{q,oj) are FDR matrix propagators then 

Dafiiq, w) = 51 ^"^(9' oj)(^fiu{q)Cui3{q, oj) (229) 

fit/ 

is also a FDR matrix propagator. The proof is rather direct. Look first at the response 
channel: 

DBB{q,uj) = ^AB^,{q,uj)aBu{q)Cufi{q,(^) 

= ABpiq, u)app{q)CpB{q,u;) = ■ (230) 

Consider next the off -diagonal componets 

DpB{q,uj) = ApB{q,uj)aBp{q)CpB{q,uj) 

= ApB{q,u)V{q)CpB{q,u) (231) 

Dspiq, uj) = Aspiq, uj)V{q)CBp{q, oo) (232) 
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It is easy to see that 

DpB{q,oo) = D%^{q,u) . (233) 

Next consider the diagonal component 

Dppiq, uj) = App{q, uj)V{q)CBp{q, uj) + ApB{q, uj)V{q)Cpp{q, uj) 
V{q) 



if3ijj 

V{q) 



[{ABp{q,uj) - ApBiq,uj))CBpiq,(^) + v4pij(g, w)(Cijp(g, w) - CpB(g, w))] 

[{ABpiq, uj)CBpiq, ^) - Apsiq, co)CpBiq, ^))] 

[DBp{q,co)-DpB{q,oj)] (234) 



if3u 

1 



if3u! 



or 



i(3ujDpp{q, u) = Dspiq, co) - DpB^q, w) (235) 

and 

- ^/3w^pp(?, w) = DpB^q, uj) - DBpiq, uj) (236) 

Together Eqs. (l235p and (I236P give 

Dppiq,uj) = D;piq,uj) (237) 

and 

D^fs{q,uj) = D;^{q,uj) . (238) 

Then 

Gafiiq, uj) = Gap{q, uj)af,u{q)Gyii{q, uj) (239) 

is a FDRMP, as is 

Gai3{q,uj) = Gap{q,u})crp^{q)g^fj{q,u) . (240) 

In operator notation 

G = gaGag . (241) 
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